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We develop the formalizm for calculating semi inclusive observables at high energy in the 
JIMWLK/KLWMIJ approach. This approach is valid for scattering of a small perturbative pro- 
jectile on a large dense target. We consider several examples including diffractive processes, elastic 
and inclusive over the target degrees of freedom, scattering with fixed total transverse momentum 
transfer and inclusive gluon production. We discuss evolution of these observables with respect to 
. . . various rapidity variables involved in their definitions (total rapidity, rapidity gap, width of diffrac- 

^>0 ' tive interval etc.). We also discuss the dipole model limit of these observables and derive closed 

' simple (as opposed to functional) differential equations in this approximation. We point out that 

, there are some interesting differences between the full JIMWLK/KLWMIJ evolution and the dipole 

model evolution of diffractive cross section. In particular we show that in the dipole approximation 
^ ' the target does not diffract beyond the valence rapidity interval, consistently with the intuition 

Q , suggested by the Pomeron fan diagramms. On the other hand such diffractive processes do exist in 

^ . the fuU JIMWLK/KLWMIJ approach, albeit suppressed by the factor 

in 

(N 



X 



I. INTRODUCTION AND CONCLUSION 



(N 
> 

' Attempts at understanding of hie;h energy hadronic scattering have been at the forefront of the QCD research for 

V/^ , a very long_time [3, H, H, H, H, [3, H, H, iU ■ In recent years these have increasingly focused on the physics of 

Ps| saturation J3| , , [llj in the formal framework of what is being referred to intermittently as the Balitsky-Kovchegov 



OO (BK) [12] , [Q or the JIMWLK evolution equations Ij] , [la] ■ This set of equations describes evolution of the scattering 



' amplitude of an arbitrary perturbative hadronic projectile on a dense hadronic target. The evolution of the projectile 



[ wave function in this approach is perturbative, but the interaction with the target is resummed to all orders in the 

■ target gluon field strength in the eikonal approximation. Recently a lot of thought has been devoted to finding the 
' extension of this approach to include nonlinear effects also in the projectile wave function - the so called Pomeron 

loops [H, 113, Hi, [3, is mi 111 H, m, m, M M, MMMMMMMMM- n is also worth mentioning 
O ^ that the JIMWLK/KLWMIJ equation and its generalization including the Pomeron loops has been interpreted as the 
D , perturbative QCD incarnation of the Reggeon field theory [sj, and some properties of its Hilbert space and spectrum 
• have been studied ^36]. 

■ The BK/ JIMWLK evolution equation directly describes the evolution of the total cross section. Indeed most of the 
work in this framework has concentrated on understanding of this quantity. On the other hand, as is now clear, the 
total cross section at asymptotically high energy cannot be reliably calculated in this approach, since it violates the 

' Froissart boundjsi]. On the other hand one does expect that other more exclusive observables have a better chance to 
be calculable in the perturbative saturation approach. It is thus important to understand how to adapt the formalism 
to calculation of such semi inclusive quantities. This is the main purpose of the present paper. We hasten to note 
that several semi inclusive observables have been already discussed in this or similar frameworks. In particular in the 
dipole model approach diffractive dissociation cross section was considered in [3^ , single gluon inclusive spectrum was 
calculated in [io] and double gluon emission spectrum was discussed in [4l|. Within the JIMWLK approach proper, 
the single gluon spectrum and the double gluon spectrum with rapidities of the two gluons close to each other has 
been calculated in [i^] and the formalism for calculating diffractive cross sections has been put forward in [4^] . 

The aim of the present paper is to consistently and fully develop the formalism for calculating semi-inclusive observ- 
ables within the JIMWLK/KLWMIJ framework. Many elements of the developments in this paper are independent of 
the explicit form of the JIMWLK Hamiltonian and those generalize straightforwardly to the as yet unknown complete 
high energy evolution Hamiltonian including Pomeron loops. However by using the JIMWLK Hamiltonian for the 
evolution of the target wave function we restrict the applicability of the current calculation to large target and small 
projectile. More specifically, the target already at zero rapidity is assumed to be a dense system, such as nucleus. 
This allows one to use JIMWLK evolution for the target wave function already at zero rapidity. The projectile on the 
other hand is assumed to be a dilute system not only at the initial rapidity but all the way to the total rapidity of 
the process Y. This allows one to use KLWMIJ evolution for the projectile wave function all the way up to the final 
rapidity Y. These conditions ensure that the "pomeron splitting" contributions to the evolution of the target wave 
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function and the "pomeron merging" contributions to the evolution of the projectile wave function are unimportant 
throughout the whole rapidity range considered. The calculation is not appropriate to "proton-proton" scattering at 
asymptotically high energies, but rather is pertinent to "DIS on a nucleus" at preasymptotic energies. We stress that 
this limitation is not specific for application of JIMWLK evolution to semiinclusive observables, but is inherent in the 
JIMWLK approach as a whole. 

Our starting point is the energy evolution of the light cone wave function [2^, [3] . We use essentially the same 
approach as in [42l | and develop it further to include a wider class of observables. We consider several examples, in 
particular variety of elastic and/or diffractive cross sections, cross section with fixed transverse momentum transfer 
and inclusive gluon spectrum. In all these cases we define the appropriate observable, derive its evolution with rapidity 
(total rapidity of the process and/or width of the rapidity gap and/or width of the diffractive interval) and discuss in 
detail the dipole model limit for each one of the observables. We make an explicit link with the work [43| and provide 
explicit definitions of the basic quantities used in [4^ in terms of physical S matrix elements. 

Apart from setting the unified framework for discussion of semi inclusive observables, we make several specific 
points which are worth noting. First, in addition to the diffractive cross section discussed in [s^ which is elastic 
in the target degrees freedom, we consider N£ which is inclusive over the final states of the target within a small 
target side rapidity interval. We confirm that the two quantities are different, and derive the dipole limit evolution 
for both (equation for being the Kovchegov-Levin equation [1^). We also show that the more general double 
diffractive cross section in the dipole approximation does not depend on the width of the diffractive interval on the 
target side, but only on the total distance in rapidity between the target and the diffractive remnants of the projectile. 
This property is straightforward in terms of the pomeron fan diagrams, since once the final state cut enters the gap, 
it cannot cross any more Pomerons, and thus the width of the diffractive interval on the target side is always confined 
to a finite fixed value. Interestingly, we show that in the full JIMWLK/KLWMIJ framework this property does not 
hold, and diffraction of the target is possible. Thus the 1/Nc corrections to the dipole model which are present in the 
JIMWLK equation lead to a qualitatively different behavior of the diffractive cross section. 

Any experimental observation of double diffractive processes can be viewed as a measure of physics beyond the fan 
diagram approximation. This physics is traditionally associated with effects of pomeron loops, which have double 
diffractive processes as experimental signature. We observe, however, that the JIMWLK/KLWMIJ framework beyond 
the dipole model limit leads to similar signatures even though it does not take into account high density effects in the 
projectile wave function. 

This paper is structured as follows. In Section 2 we recall the general framework of the light cone wave function 
evolution and how it leads to the JIMWLK/KLWMIJ evolution equation. We also discuss the large Nc limit and how 
it leads to the dipole model and its generalizations. In Section 3 we introduce the method of calculating semi inclusive 
observables in this framework and consider variety of elastic and diffractive observables and derive their evolution 
equations with respect to total rapidity as well as rapidity gap and diffractive rapidity intervals. In Section 4 we derive 
the dipole limit of the evolution equations in all the cases. This is necessary for any possible numerical implementation 
of the evolution with the view of phcnomcnological applications since it reduces the evolution equations to simple 
rather than functional differential equations. In Sec. 5 and 6 we discuss the scattering cross section with total 
transverse momentum transfer and the inclusive single gluon production, as examples of non-diffractive observables 
which are calculated by the same methods. 



II. HIGH ENERGY SCATTERING: THE GENERAL SETUP 

A. Evolution of the wavefunction 

In this section we review the approach to high energy evolution based directly on the evolution of a hadronic wave 
function. We will follow the formalism of [i^. We concentrate on the gluonic part of the wave function, although 
including quarks does not pose any additional differences. 

Consider an arbitrary high energy hadron with large rapidity Yq. Its wave function in the gluon Fock space can be 
written as 

\v)yo = *k:'.(^)]|0>yo (II.l) 

The gluon creation operators depend on the transverse coordinate, and also on the longitudinal momentum 
The gluon operators in eq. pi.l|) all have longitudinal momenta above some cutoff A. We refer to these degrees 
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of freedom as "valence". Henceforth we omit the dependence on longitudinal momentum in our expressions, as the 
momentum enters only as a spectator variable and only determines the total phase space available for the evolution. 

When boosted, the valence wave function gets dressed by the evolution "Cloud" operator fl Physically this 
operator creates the soft Weiszacker- Williams gluon field which accompanies the boosted valence gluons. It therefore 
depends on the color charge density of the valence degrees of freedom (p) and also involves creation operators of softer 
modes ak+ with longitudinal momenta below A. The evolved wave function has the following structure 



I*) 



nvip, a) \v)yo 



(11.2) 



The evolution operator is explicitly known in the dilute limit (p 0) only. In this limit, for Y — Yq = SY <C 1, it 
is the coherent operator 



^p^o = CsY = Exp <' i / (fxb1{x) 



A 7ri/2|fc+|i/2 



(11.3) 



Here the creation operators a\k'^) create gluons with soft momenta, which are not present in the valence state \v) 
The field b depends only on the valence degrees of freedom through 



2-K 



, [z - X)i 
[z — xY 



p^{x) . 



For a finite evolution interval the appropriate expression is 



-PExp 



d^xbf(x,k+) — 



dk+ 



7rl/2|j!c+|l/2 



<(fc+,x) + al''{k+,x)) 



(11.4) 



(11.5) 



where V denotes the path ordering with respect to k'^ , and the k~^ dependent field b includes the charge density of 
all modes harder than fc"*" 
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, {z - x)i 
{z — xY 



dp+a^(p+)T"a(p+)] 



(II.6) 



One can resum some corrections away from the low charge density limit by keeping the full nonlinear relation between 
the Weizsacker- Williams field and the color charge density [1^ |3l[. In this case b^ is determined as the solution of the 
" classical" equation of motion 

d,b1 + ge'^''b\{x)b1{x) = gp%x) 

e^[d,b'; - d,bt + ge^'^lb^^] = (H.7) 
The coherent operator C dresses the valence wave function by the cloud of the Weizsacker- Williams gluons: 



C^At{k+,x)C = At{k+,x) + —bnx). 



Given this evolution of the hadronic wave function one can calculate the evolution of an arbitrary observable 0{p) 
which depends on the color charge density. For example, the S matrix in eikonal approximation belongs to this class 
of observables. In the by now standard notation 

{v\6[p]\v) = I DpW[p] 0[p]. (II.9) 

As discussed in detail in [l^l, the integration variable p on the right hand side of this equation is understood to 
depend on transverse coordinates as well as an additional coordinate a;~, which can be thought of either as the 
longitudinal coordinate of the hadronic wave function or as a mathematical "ordering coordinate" present to enforce 
correct commutation relations of the operators p'*. 

The evolution of the expectation value is then given by 

d{v\d\v) {v\n^yd{p + 5p)nY\v) - {v\d{p)\v) 

-^Y- - y'^^Yo Y-Yo 



J DpW[p]H[p]0[p]. (11.10) 
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The color charge density in the first term contains also the contribution of the soft gluons generated by the evolution 

.e^-^o A 

Sp^ix) = / dk+al\k+,x)T^,a'i{k+,x) (11.11) 

J A 

where — if'^^'^ is the SU{N) generator in the adjoint representation. This is conveniently represented in terms of 
the charge density shift operator (which also has the meaning of the "dual" to the Wilson line operator) 



■ 1 



R{z 



The evolution Hamiltonian H generally can be written as 



Sp^iz, z ) 



(11.12) 



S d 

HiP,^] = - ^ {Oa\nl{p,a)nYiRp,Ra)\Oa}\Y=Yo = Qn[R] Q n[R] ■ (H-IS) 

n 

The last equality in pi.l3[) is given in terms of the n-particle production amplitudes Q„ [3^. Each Q„ depends on 
n transverse coordinates, color and Lorentz indices which we do not indicate explicitly. This expression is formally 
valid for arbitrary charge density. To write down an explicit expression for H we specify to the dilute regime. In this 
case only one soft gluon is created in one step of the evolution and the Hamiltonian gets contribution only from the 
one gluon (n — 1) production amplitude: 

HKLWMij^pj_^ ^ ^ A. (^0^\Cl{p,a)CY{Rp,Ra)\Oa)\Y=Y„ - I ^Q1{z)Q-{z) (11.14) 

where the one gluon emission amplitude Qfiz) is defined as (we suppress the index n — 1) 

Qf{z) = R'"'{z)b1,{z) - hl,{z) (n.l5) 
Here the gluon field b is the function of the SU{Nc) rotation generators 

KLim(^^) ^'-{^ -JUmi-) ■ (11-16) 



with 



r^(x)=-tT\Rix)T^jj^Y j£(2,.)^_tr|™(x)^^^, Jlix) = [R{x)Mxr. (H.17) 



so that 

^'^^'^ ^£ J '^'''^^ [^"'^^^ ~ ^"'("^^i -^^^"^^ ■ ^"-^^^ 

We refer to the evolution Hamiltonian in this dilute limit as the KLWMIJ Hamiltonian ^2^. 

The partial resummation of the nonlinearities mentioned above which corresponds to keeping 6L.fi as solutions of 
the full classical equations of motion with sources given by the generators of the left/right color rotations Jr. B[25(.i3Tj| 
that leads to the evolution Hamiltonian hI<lwmij+^ 

We note that for the derivation of eq. (jH.lOp it was not crucial to consider a diagonal matrix element of the operator 
O. The same derivation can be repeated straightforwardly for a generic non-diagonal matrix element as well. Defining 

{v\d[p]\v') = J DpW[p] 0[p] (H.19) 

we find the evolution 

d{v\d\v') ,. {V\n\.d{p+5p)ny\v') - {v\d{p)W) / ^rr,^ n( M '\ 

DpW[p]H[p] 0[p] (H.20) 



with the same Hamiltonian as in ea. pi.l3|l . We will need to use this fact in the following. 
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B. High energy scattering 



Throughout this paper we treat the scattering of fast gluons of the projectile on the target in eikonal approximation. 
We denote a S-matrix of a single gluon which scatters on a fixed configuration of chromoelectric field of the target by 

S^''{x)^{0\atix)Sc4\xm (11.21) 

where S is the second quantized S'-matrix operator of the field theory which in the eikonal approximation [g^l is 



exp<i / d xpp{x)a^{x) 



(11.22) 



In the natural projectile light cone gauge {A~ — 0) the large target field component is , which we denote by ax- 
The single gluon S-matrix S"''' does not depend on the polarization of the gluon and is diagonal in the transverse 
coordinate x and is given by 



S{x) = V exp{i [ dy-T''ar{x,y^)}. 
Ja 



(11.23) 



For a composite projectile which has some distribution of gluons in its wave function the eikonal S'-matrix can be 
written in the form analogous to S{x), see [1^ 



(P\S\P) 



dppW [pp] expii dy dxpp{x,y )a^{x,y ) 



(11.24) 



with Xi - transverse coordinate. The operator pp{xi) is the color charge density in the projectile wave function at a 
given transverse position, while T/F^[p] is the same weight functional as appearing in eq. (|ILTo1) . For a single gluon 
p^{xi) = T''5'^{x^ - x°), and eq. ((TL24)) reduces to eg. ([1123)) . 

To obtain the total S'-matrix of the scattering process at a given rapidity Y one has to average E of eq. (jll.24l) over 
the distribution of the color fields in the target. Denoting the probability density to find a particular configuration 
aT{x,x~) by W'^[aT{x,x~)] we write 



S[Y) 



Da^ Wl[aT{x,x-)] ^(.^YoMx^x-)]. 



(11.25) 



In Eq. pi.25p we have chosen the frame where the target has rapidity Yo while the projectile carries the rest of the 
total rapidity Y — Yq. Lorentz invariance requires S to be independent of Yq. 

The high energy evolution of the S-matrix follows from cqs. pi.lOIII.24p as 



= - j Da?r W^^[aT{x,x-)] H 



S 



5 aT 



^Y-Yoio^Tix^X )] 



(11.26) 



with the Hamitonian eq. (jII.13p with p substituted by ^^.a^'f ^-^ The Hamiltonian can be viewed as acting either to 
the right on E or to the left (as it is Hermitian) on W: 



_d_ 



^pp 



H 
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^pp 



_d_ 
W 



H 



ax, 



5ai 



W'^[aT] 



(11.27) 



As was shown in |26| in order for the total S-matrix to be Lorentz invariant and symmetric between the projectile and 
the target, the evolution Hamiltonian H must be self dual[l^. That is it has to be invariant under the Dense-Dilute 
Duality transformation 



Q°'{x,X ) 



Sp"'{x,x ) 



da'^{x, X 



-ip^-ix.x ), 



S ->R 



Hence 



H 



5a 



Y^QiiS] Qn[S]. 



(H.28) 



(n.29) 
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In the situation where the target is large and the projectile is small the symmetry between the target and the 
projectile is irrelevant. In this limit the Hamiltonian is given by the JIMWLK expression [3, [IB| which is the dual 
of the i/™M7j 



^JIMWLK ^ / Qa(^^)Qa^^) (11.30) 
J z 

with Q[S] obtained from Q[I{\ (jII.15[) by substitution S for R. 

Some of the derivations in this paper are independent of the explicit form of the Hamiltonian. However whenever 
the explicit form is required we are going to use the JIMWLK Hamiltonian cq. (jll.30p with 

Q^(z, [S]) = ^ J d'x^^-^ [S^'iz] - S^'ix)] Jlix) (H.31) 

where 

J^x) = -tr !^S{x)T<^j^^ , Jlix) = -tr |r-5(a:)^^| , J^ix) = [S{x) Mx^. (H.32) 

We record here two properties of this Hamiltonian which will be useful in our discussion of the evolution of diffractive 
observables. For any function F 

[Qtiz,[S]) + Qtiz,[S])]' FISS^] = [ [Q^zASS])]' ^SS] 



[QUz, [S]) + Qtiz, [S])y F[S, S]\ = [Qtiz, [S])f {F[S, S]} . (IL33) 

s=s 



where S and S are both arbitrary unitary matrices. 



C. Large Nc lore: dipoles, quadrupoles and such. 



The general setup for the high energy evolution is that of functional evolution equations for the scattering ampli- 
tudes, or equivalently effective quantum field theory for the unitary matrix S. Some aspects of this Reggeon field 
theory have been studied in 'st'Ij^g']. The problem of its solution is however a formidable one, even though the theory 
is considerably simple than full QCD. It is thus desirable to have a simple truncated version of the theory which would 
reduce the complexity below the level of quantum field theory. Such a truncation is offered by the formal large Nc 
limit and its simplest variant is Mueller's dipole model [gI. |47|. 

Any physical projectile must be color singlet, and the simplest color singlet state is a fundamental dipole. Assume 
for a moment that the projectile wave function contains only dipoles as in Mueller's dipole modelQ. How does the 
scattering matrix of such a projectile evolve? In the large Nc limit color singlet objects evolve independently of each 
other. Thus every dipole in the wave function leads his independent life. Formally this means the following. The 
scattering matrix of a single dipole is 

six,y)^^tr[SFix)sUy)] (11.34) 

where the subscript F denotes fundamental representation. The S matrix of a dipole projectile is therefore some 
function of the variable s only 

E^^[5] =S^^[s]. (H.35) 
For a given projectile wave function the S'-matrix S[s] is easily calculated as 

E^^[s] = ^ P„{(xi, yi), {xn,yn)}s{xi,yi)...s{xn, Vn) (11.36) 

n 

where P„{(xi, yi), (a;n,yn)} is the probability to find n dipoles at the specified points in the incoming projectile 
wave function. 
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The high energy evolution of such a wave function in the large Nc limit obeys Mueller's dipole evolution re-expressed 
as Hamiltonian evolution in The same result is derived starting directly from the JIMWLK equation [25]. The 
result is 

dY 



with 



with the dipole kernel 



^dipolc = |1 /■ M,,y,,[s{x,y)-six,z)siy,z)]—^; (11.38) 



- yf 
{x — zY (z — yY 



M.,y,. ^ ,^ (11.39) 



The solution of the dipole evolution equation eq. Eq. pi.37p can be expressed in terms of the solution of a simple 
differential equation [48| 

SF[s] = ^?-nsY] (11.40) 



where sy solves the BK equation (see Refs. [50|, |51|, |52|, |53|, |5J| for analytical and numerical studies of the BK equation) 

Mx,y^^[sY{x,z)sY{z,y) - SY{x,y)] (11.41) 



dsy 



with the initial condition 



SYo = s (11.42) 
At the same time the expression for the probability density, Eq. (|II.25|) , turns into 

S'''^°'^{Y) = I Ds <H E^^yJ.], (11.43) 

i.e. the average over the gluon field of the target is rendered as an average over an ensemble of dipoles at the 
initial Yq used in Eq. pi.4ip . This average still allows to accommodate nontrivial, non-factorized multi-s correlators 
{s{xi,yi)---s{xn,yn))T, see Hi in Si [^i^. 

Further simplification is achieved if one assumes that the dipoles scatter on the target independently. This amounts 
to factorization of the target averages of the dipole s-matrices 

{s{x,y)s{u,v))T = {s{x,y))T{s{u,v))T (11.44) 

With this assumption, one replaces the ensemble average over target fields, or alternatively over the ensemble of 
functions SYoix,y) shown in pi.43[) . with a fixed initial function SYoix,y). We will refer to this factorization property 
as the target mean field approximation. Within the target mean field approximation 

(^rMr = sr[(^r>T] . (11.45) 

We stress that this mean field approximation does not follow logically from the dipole model approximation for the 
evolution kernel eq. pl.37p . but rather is an additional assumption about the properties of the target. Physically it 
means that the target fields are assumed to fluctuate so strongly that they are completely uncorrelated in different 
points in space. In practice this is good assumption if the points are separated by a distance larger than the saturation 
length Rs = Qj^, which is also the correlation length of the target fields. However for two dipoles separated by a 
distance smaller than Rg in the impact parameter space this approximation is bound to break down. We will come 
back to this point later. For a more detailed discussion see [Zsj . 

The dipole model provides the simplest known framework for model discussion of high energy evolution. However 
as we shall see below, for some observables it is not sufficient to consider dipole degrees of freedom alone. In this 
case we have to allow for existence of quadrupoles. Fortunately the large Nc approach generalizes effortlessly beyond 
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the dipole model. It is not necessary to assume that the projectile contains only dipoles. One can indeed allow the 
function to depend also on the quadrupole degree of freedom 

q{x,y,u,v) = ■^tr[SF{x)Sl{y)SF{u)Sl{v)]. (11.46) 

The crucial property of the large Nc evolution is that all singlets evolve independently. It is then straightforward to 
show that the evolution of a function Yj^^[s. q\ in the large Nc limit is 

'^^y']^^'-''^^ = _ ( ijdipolc ^ ^quadrupole) ^PP ^] (jj 47) 

with ijd^o'e given by eq. (|TL37)) and 

^quadrupolo _ ^-jj ^g) 



27r 



x,y.u,v,z 



Lx.v,u,y;z ^x,y ^u,v H~ L/x^v ,u,v;z Qx,y,u,z ^z.v ~^ -^x.y ,x ,v.z Qz^y.u,v ^x^z 
^~ ^x^y ,u.y\z Qx,z,u,v ^z^y H" -^u^y^u^v^z Qx^y,z,v ^u,z} 7 7 (11.49) 

oq[X^ V) 



where 



-^x,y,u,v;z 



{x - z)i [y - z)i 



[x — z)^ {y — zY 



{u - z)i (v - z)i 



(u — z)^ (v — zY 



(11.50) 



Again, the solution for the evolution equation ca. pi.47[) reduces to solution of ordinary (as opposed to fmictional) 
equations 

E^ts,?] =sr[5r,gy] (11.51) 

where the dipole sy solves the Kovchegov equation eQ. pL4ip and the quadrupole solves an analogous quadrupole 
evolution equation (see Appendix A and Ref. f4ll |) 



dQx,y.u,v 



{ [^'^x,y;z ^u,v;z -^x^u,v,y;z] Qx,y.u,v -^x.,y,u.v:z ^x,v '^y,u -^x,v,u,y;z ^x.y ■^u^v~\~ 



dY 2tt 

~^ .v,u,v;z Qx.y,u.z ^z,v ~i~ -^x .y ,x ,v\z Qz,y,u,v ^x,z ~l~ -^x.y,u,y;z Qx,z,u,v ^z.y H~ ^u,y ,u.v:z Qx.y,z,v ^u,z} 

(11.52) 

with the initial conditions 

SY=Q = s] qY=o = q- (11.53) 

In principle for this mixed dipole-quadrupole model one can again apply the target mean field approximation by 
assuming factorization of all the singlet averages 

{s{x,y)s{u,v))T = {s{x,y))T {s{u,v))t; (11.54) 
{q{x,y,u,v)s{z,z))T = {q{x,y,u,v))T {s{z, z))t ; 
{q{x,y,u,v)q{x,y,u,v))T = {q{x,y,u,v))T {q{x,y,u,v))T ■ 

Generalization to higher multipoles is in principle straightforward, but we will not need for the observables considered 
in this paper. 



III. SEMI-INCLUSIVE REACTIONS 
A. Generalities 

We are interested in calculating characteristics of the final states emerging after a collision of the projectile, which 
at the initial rapidity has the wave function and the target with the wave function \Ty). The target and the 
projectile are boosted before the collision to the total rapidity Y. 
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Let ily^Yo ^'^^ "i^n-ots evolution operators for projectile and target boosted to the rapidity Y — Yq and Yq 
respectively. The time of interaction is set to be i = 0. 

The total wave function coming into the collision region is therefore at time t = 

- nP_y^ |p„) |r„) . (iii.i) 

The system emerges from the collision region with the wave function 

|*out) = sn^_Y, n^, \Pv) \n) . (111.2) 

The system keeps evolving after the collision to the asymptotic time t — > +00, at which point the measurement of an 
observable O is made. As explained in the evolution of the outgoing wavefunction from t = to t — > +00 is given 
by the Hermitian conjugate of the same operator f2. Thus the general setup for computing any observable O in the 
final state is 

id) = {n\ { P.I np^y^ f]^; (1 - ^t) qp_^^ f^pt t _ ^P_^^ (jjj_3) 

where the factor S — 1 ensures the proper subtraction of the contribution of the initial state. Generically, the 
observable O may depend both on target and projectile valence degrees of freedom as well as gluon degrees of 
freedom at intermediate rapidities. An example of rapidity dependent observables is diffraction considered in the next 
subsection (for example eq. (|III.29p and single inclusive gluon production discussed in Sect 6. 

To express this in terms of scattering amplitudes we insert the resolution of identity on the projectile and target 
Hilbert spaces 1 — \p') {p'\ and 1 — \t') {t'\. Then we have 

(6) = {TyJS^P' - \t') 01,^; {t"\5^"P - ti:\ \Ty,) (III.4) 

where we introduce a non-diagonal S-matrix element which is an operator acting on the target's Hilbert space 

S^!_Vo = {Py-y^S\p'); (III.5) 

and the matrix element for the observable 

of,;^: ^ {t'\ {p'\n^_y^ nl 6 n^yly^ nl^ \p") \t") . (iii.e) 

It is important to keep in mind that the resolution of identity on the projectile side is achieved by full basis of 
states in the Hilbert space of the projectile and not only by those states that can be obtained from valence states 
evolved to the projectile rapidity. In this sense the incoming states are of a very special type, as by definition we only 
consider evolved valence states. The intermediate states for example include states which have no valence gluons at 
all, but contain only soft gluons. Such states do not appear as initial states for the processes we consider. To make a 
distinction between the incoming states and the basis states that span the full Hilbert space, we denote the former by 
capital letters while the latter by script letters. Thus for example the ^-matrix element between two evolved states 
will be denoted by 

E^^;„ ^ (P.|f^^ly„ S f^^_yJP^) (III.7) 

omitting the rapidity index whenever we feel it should not cause confusion. 

A similar representation can be developed for the target side, see Appendix B. However we will not consider the 
most general observables on the target side. Instead we will concentrate on processes which are either completely 
inclusive over the target degrees of freedom, or in which the target scatters elastically. For the processes inclusive 
over the target, the observable O does not depend on the target degrees of freedom. Thus 

The external average over target evolved to rapidity Yb can be done using the target weight functional Wy^ evolved 
from zero rapidity with the Hamiltonian jjJIMWLK ^ 

(o)ti = j Dswlis] J2 {s'''' - ^l-%,is]) {s""^" - s^p^Vj^]) o^'p". (in.9) 

p'p" 
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For the target elastic processes only one intermediate state, namely \T) contributes. In this case we get 



(O) 



TE 



(Ty„\ [S^P' - ±1-%^) \Tr,) O^'P" (TyJ [S^P" 



y.p"P 



ITYo] 



(III.IO) 



This expression involves two independent target averages and can be written in terms of the target weight functional 
W as 



(O) 



TE 



SPP" 



qp: 



(III.ll) 



In these equations we have written the integration measure as DS rather than Da as in eqs. (jII.25|II.26p . This we do 
purely for notational simplicity. The actual integration measure is indeed Da, but since the integrands that we will 
encounter depend only on S we allow ourselves this notational shortcut. With this note of caution we will use this 
notation throughout the paper. 

Within the above formulation, the problem of high energy evolution of any observable reduces to identification of 
corresponding quantum operator O and its matrix element C ^ . 

In the rest of this section we will discuss several observables related to elastic and diffractive scattering and their 
evolution with respect to various rapidity variables on which they depend. To make clear our notations, superscript 
over any observable refers to the projectile and subscript to the target. Thus for example N^'^^ denotes the cross 
section for projectile diffraction with the diffractive interval Yp with target scattered elastically. The total rapidity 
of all processes discussed here is denoted by Y and most of the time we will not indicate it explicitly. 

We start our discussion with simple observables which depend only on total rapidity. 



B. Elastic Scattering 



P P 



P P 



T I T 
(a) 




(b) 




P' P 




P' P^ 



T T 

(c) 



(d) 



FIG. 1: Elastic scattering: (a) total elastic: the final states on the cut are the same as the initial ones both on the projectile 
(P) and target (T) sides; (b)projectile elastic, target inclusive - all final states on the target side T' are summed over; (c) target 
eslastic, projectile inclusive: all projectile final states P' are summed over; (d) double inclusive with maximal gap. Horizontal 
lines denote final state gluons. 



1. Total elastic scattering 

Consider a totally elastic process where both the projectile and the target scatter elastically Fig. [Ua. For the target 
averaging we use eQ. (jIII.ll[) . On the projectile side the observable O is 



The total elastic cross section reads: 



o 



\Pv) {Pv\ 



I DSWliS] (1 - E^^^J^]) 



(III.12) 



(III.13) 



This expression obviously does not depend on the rapidity Yq which separates the projectile and target Hilbert spaces. 
The evolution of the amplitude Ty = 1 — Sy is given by the JIMWLK equation eq. (|lL26)) . 
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2. Projectile elastic scattering inclusive over the target 



This observable (Fig. (Ub) is defined as elastic over the valence degrees of freedom of the projectile and inclusive 
over all the rest of the rapidity interval. Thus we should take the separation rapidity at the total rapidity of the 
process Yq — Y. Using ea. pil.9|) with the operator O given by eq. piI.12[) we get 

Nf ^ J DSW^[S] |l - ^o^[S]f . (III. 14) 

The evolution of this observable is given by the JIMWLK evolution of Wy [S] . 



3. Elastic target inclusive over the projectile 



Target valence degrees of freedom scatter elastically while all states on the projectile side over the full rapidity 
interval are summed over inclusively (Fig. [He). We use eq. (|III.llj) 

Ne = J DSDSW;i^[S]W^[S] (1 - - ^^^[S^ + (III.15) 

We remind the reader that we are interested in the situation where the projectile is dilute and the target is dense 
all the way through the evolution. This allows us to choose Yq = in ea. piI.lSp . The evolution of each S term in 
eg. dlll.lSp is given by the JIMWLK Hamiltonian which depends on the argument of S. The term of most interest is 
the last term involving Sy^ [S*^ S] . Using the property of the JIMWLK Hamiltonian eq. (jll.33|) we can be written 

^ s] = -H3[S; S] S] (III. 16) 

We have introduced the operator [11] 

Hs[S;S] = / [Qn^dS]) + Q-{z,[S]t ■ (in.i7) 

J z 

The Hamiltonian in this equation can then be interpreted as acting to the left on the target weight functionals. Thus 
the amplitude can be represented in the following form 

Ne = J DSDSZy{S,S) (1 - - S^^[5t] + ^^^[3^3]) (IH.18) 

with the functional Z satisfying the functional equation 

^Zy[S, S]^ - HslS; 3] Zy[3, S] . (HI. 19) 

The initial condition for this evolution is 

Zy^o [S, S] = [3] [3] . (III.20) 



C. Diffraction 



An interesting set of observables is that of diffractive observables with rapidity gap. We will be interested in 
identifying the evolution governing the diffractive observables both with respect to the total rapidity of the process 
Y as well as with respect to the rapidity gap(s). 
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FIG. 2: Diffraction: (a) Projectile diffraction; (b) Target diffraction; (c) Double diffraction; (d) Central diffraction. Notations 
are the same as in Fig.l. 





P P P P 





(b) 



(c) 



FIG. 3: Fan diagrams for diffraction: (a) Projectile diffraction with target scattered elastically; (b) Projectile diffraction with 
target diffracting in a small rapidity interval; (c) Projectile scatters elastically. Notations are the same as in Fig.l. 



1. Double inclusive with maximal rapidity gap 



The simplest observable of this type is the cross section inclusive over the valence degrees of freedom of both 
the target and the projectile and with rapidity gap covering the whole rapidity interval Y (Fig. [l]d). We put the 
separation rapidity Iq close to the target, Iq = 0. From the target point of view this is an observable of the type 
eq. (|IIL9|) . From the point of view of the projectile we have to understand which intermediate states are allowed to 
contribute in the sum over p' . Since we are requiring that no soft gluons be found in the final state, clearly the 
only states that can contribute are the states of the type P, that is projectile states which can be obtained from a 
valence state by boost to rapidity Y . Moreover all such states should be summed over with equal weights. Thus this 
observable is given by 

= / DSW^[S] Y: ('5^^' - ^V^'iS]) {S""" Sr^[5]) . (111.21) 
J p, 

The evolution is simply given by the evolution of each factor 'S[S] according to the JIMWLK equation. This again 
can be recast in the form in which the gap as well as the evolution is attributed to the target degrees of freedom. To 
do this we rewrite the evolution for eq. (jIII.2ip as follows 

-^K:o = - I DSDSW^iS] S{S-S] X (III.22) 
p' 

with 

H2 [S; S] = i/^^A^M/LX j^j ^ jjJIMWLK j^j gjj_23) 

Now understanding the action of the evolution operators to the left rather than to the right we can write the same 
observable as 

Nofi = / DSDS (1 - Sri^] - Sri-?^] + ^o'^iS^S]) Zy[S, S] (III.24) 
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with Zy [S'j S\ satisfying the evolution equation 

^Z\SA = - {u^nMWLK^g^ ^ Z{S,S\ (III.25) 

with the initial condition 

Zy=q = W'^\S\ 8\S~'S\. (III. 26) 

To write ea. pil.24l) we have used the fact that when the projectile wave function is un-evolved, the states P' constitute 
complete basis in the projectile Hilbert space, and that the projectile averaged S- matrix E considered as a matrix 
on the projectile Hilbert space has the property 

= YPP\S''S\. (III.27) 

We now consider processes where the projectile diffracts into a rapidity interval Yp . This interval is not necessarily 
small, so this type of observable can be evolved independently over the total rapidity Y and the width of the diffractive 
interval Yp. The target can either scatter elastically or can in principle also diffract. We start with the process where 
the scattering on the target side is elastic, the process considered in (39j . 



2. Projectile diffraction with target scattered elastically 



Consider the projectile diffraction with rapidity gap imposed on the target side where the target undergoes elastic 
scattering (Figs, dja and[31a). The diffractive interval on the projectile side is Yp and the rapidity gap Ygap = Y — Yp. 
We choose the rapidity separation Yq to be at the target rapidity Yq — 0. We show in Appendix B that the result does 
not depend on the position of the separation rapidity, as it should. Clearly we should use eq. pil.lip on the target 
side. On the projectile side we have to sum over all intermediate states which can be obtained by the evolution over 
the gap. Thus the intermediate states p' in eq. (jIII.ll[) should have the form 

\p) - ny,^, \PY,) (111.28) 

where the states \pYp) form a complete basis inside the rapidity interval Yp and not just on the valence part of the 
projectile Hilbert space 

In the framework of eq. (|III.3|) this corresponds to choosing the operator O = Sp^^ with \X) — fly^lpYp). 

Note that the states \X) defined through this relation are not necessarily valence states, as remarked in the discussion 
in Section HI. A. In fact, since \pyp) span the full Hilbert space on the rapidity interval Yp, so do the states \X), since 
i^Yp is a unitary operator. Thus S]{x}|-'f)(-'^| ~ '^{pYp}\PYp){PYp\ — ^Yp where lyp is the unit operator inside the 
rapidity interval Yp. 

With this definition of the intermediate states the expression for this observable is 



D,Yp 



J DSDsw^s] w^s] (^1 - En^]- sr''[^] + E4f:[>^p;^]sCj^p;^]) • (in-29) 



Here 

^rL[^p; s] = (nif^t nl^^Jny^jPYp) ■ (in.so) 



To bring this expression to a simpler form we note that the evolution of [Yp; S] with respect to Ygap is given by 
jjjiMWLK ^pj^^g Qg^jj again integrate the evolution by parts in eq. pil.29|) and express the integrand in terms of 
T,Y^ ^q[Yp; S]. However at Ygap = 0, the states p form a complete basis in the projectile Hilbert space and we can 
use again the property eq. (|III.27p . As a result we get the target weight functionals evolved through the gap: 

TVi?'^- - J DSDS (1 - E^;[5] - Ef;[5t] + E^Jl^t^]) WlJS]WlJS]. (III.31) 
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This again can be written in terms of the target "weight function" which depends on two unitary matrices and two 
rapidities 



J DSDS (1 - E^iS] - ^riS^ + ^riSS^]) Zy,,y,JS,S] . (III.32) 



The evolution of Z with respect to the width of the diffractive interval as 

d 



dYf 



Zy,,y,., [S, S] = -H, [S; S] Zy,.y,^, [S, S] (III.33) 



with H3 [S, S] defined in eq. (jIII.17|l . The evolution with respect to the width of the gap is in general more complicated, 
however for Yp = it reduces to 

^ Zy,=o.y,^, [S, S] = - H2[S; S] Zy,=o.y,^, [S, S] (III.34) 



BY 

^ ^ gap 



Thus to find ZYp^Yg^p on^ has first to evolve Z from Ygap = to Ygap with eq. (jIII.34[) starting with the initial condition 

Zy,=o;i:,„,=o[5,5] = W^[S] W^[S] (III.35) 

and subsequently evolve the solution with respect to Yp via eq. (jIII.33p . 
This reproduces the result previously derived in Ref. [ist - 

3. Projectile diffraction with the target diffracting in a small rapidity interval. 

Another interesting diffractive observable is the cross section summed inclusively over the valence target excitations 
(Fig. Ob). Like before we fix the diffractive rapidity interval on the projectile side and the rapidity gap on the target 
side, but sum inclusively over the possible target valence states. In view of our discussion in the previous subsection 
the result is easy to write down 

= J DSW^S] (^l - St,^^[^] +5:42[rp;5]E?f^^[yp;5]j (III.36) 

This can again be rewritten in the form similar to eq. (jIII.32p 



J DSDS (1 - J:r[S] - ^ns^ + ^ns^ S]) Zy,,y,JS,S] . (III.37) 



The evolution of the functional Z with respect to the width of the diffractive interval is as before 

d 



dYp 

and its evolution at vanishing Yp with respect to the rapidity gap is 

d 



Zy,.y,JS,S] = -H3[S;S] Zy,^y,JS,S]. (III.38) 



Zy,=o,y,^, [S, S] = -H2 [S; S] Zy,=o.y,^, [S, S] . (III.39) 



The only difference relative to the observable in the previous subsection is in the initial conditions. One has to solve 
ea. pil.39l) with the initial condition 

Zy,=o;Y,^,=o[S, S] = W;[[S] 6[S - S] (III.40) 
and feed the solution as the initial condition into eq. (|III.38[) . 
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4- Double diffraction and more. 

In the previous subsection we have considered an observable which is summed over the final states of the target in 
a small (valence) rapidity interval. The restriction on the size of the diffractive interval on the target side Yp can be 
straightforwardly lifted, and we can consider double diffraction with three independent rapidity intervals Yp, Yt and 
Ygap denoted by. Note that in the framework of the Pomeron fan diagrams (Fig. [3|), this observable does not depend 
on Yt and Ygap separately but rather only on their sum Yt + Ygap. which is the distance between the valence target 
and the diffractive remnants of the projectile. This is so since the Pomerons can not be cut horizontally without 
loosing powers of energy. Thus in the absence of Pomeron loops, once a gap is required, no Pomeron can be cut 
below the gap and no particles are produced. We will show in the next section, that this property holds in the dipole 
model approximation to JIMWLK evolution. However in the full JIMWLK framework this is not the case and double 
diffraction is indeed possible. 

We start with defining the cross section where the projectile states are summed over inclusively over the valence 
rapidity, target scatters diffractively in the rapidity interval Yt and require the rapidity gap Ygap. By the same 
reasoning as in the previous section we have 

nS:1 = Jds wi [s] [i - ^z. [^1 - ^'yZ [^1 + E ^Z"^ [^] <i t^]) • ("1-41) 

Note that the sum over P' in the last term does not give unity, since P' do not constitute complete basis of states in 
the rapidity interval Ygap. 

Rewriting this in a form similar to the previous subsection we obtain 



D,Yt 



I DSDS (1 - J^riS] - ^nS^] + ^nS^S]) Zy,^^,yAS,S] ■ (111.42) 



The evolution of the functional Z with respect to the rapidity gap is given by 

d 



Zy.^^.Yr [S, S]^-H2 [S; S] Zy,^,^y^ [5, S] . (III.43) 



and its evolution with respect to the diffractive interval Yt at vanishing Ygap is 

d 



gY^ZY,^,=0,Y^ [S, 5] = - i?3 [S; S] ZY,^,=0,Yr [S, S] . (III.44) 
To calculate Zy^^^ one has to solve eq. (jIII.44[) with the initial condition(65j 

ZY,^,=o.Yr=o[S, S] = W^[S] 6[S - S] (III.45) 

and then evolve Z with respect to Ygap according to eq. (jIII.43[) . The solution can be formally written as 

Zy^.^.y^ [S, S]^e-"^ e - Zy,^,=o.y^=o[S; S] (III.46) 

We can now generalize the previous consideration to an observable with arbitrary number of gaps and diffractive 
intervals. It is clear from the discussion so far that the evolution across any gap is given by the Hamiltonian H2 
while the evolution across a diffractive interval is given by H3. Consider the process with n rapidity gaps Yi, ... , Yn 
separated by diffractive intervals yo, ... ,?/„ where yo is the diffractive interval of the target and y„ that of the 
projectile. Our previous discussion allows us to write down this observable in the following form: 



?(yi,...,r„;yo,...,2/n) = J DSDS (l - S^I^] - ^^[3^ + ^^[3^ S]) Z[Y},{yMS^ S] . (III.47) 
where 

Z[Y},{y,}[S, S] = e- e ^" ...e^^'y^ e-"'^' e'"' Zo[S; S] (III.48) 

with 

Zo[S; S] = W^[S] 6[S-S] (III.49) 
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5. Connection with the approach of fl^I . 



We briefly discuss the relation of our approach to that of [43] ■ As we have noted above, our results coincide with 
those of [43| where they overlap. To calculate the diffractive cross section with the diffractive interval Yp on the 
projectile side, the authors of [43| introduce the gluon cloud operator in the form 



C/=Ui[S,^]Uf[S,^] 



(III.50) 



with 



Ui[^,C] = Py2 exp 



dyidy2 0{yi 



(III.51) 



Uf[5, = Py, exp [i J dy^dy2 - y2)C( Jl[1])(s;;;% + . (III.52) 

Here ^ and S are both random noise variables. The correlator of ^ is identical with the the vacuum average of free 
gluon creation and annihilation operators: 

i&'y^J = {0\k+a-^\k+,u)a^''^^{p+,v)\0)=6'^''S^^6l^lSy,,y, 

The correlators of S are identical, but only defined in the rapidity interval Yp close to the projectile. 

The scattering amplitude with arbitrary number of gluons in the final state is calculated by averaging U over 
the random variable ^. The diffractive cross section is then obtained multiplying the amplitude at a given S by 
the complex conjugate amplitude at S and averaging over the random variable S, always contracting one S in the 
amplitude with one S in the conjugate amplitude. This expression is then averaged over S and S with the same 
weight functional. 

The procedure outlined above parallels exactly our approach described in this section. The equivalent of the cloud 
operator U in [i^ is the operator fly^ SQy eq. pil.3p with the operator Q defined in eq. (|II.5p . The random noise 
variables ^ and S emulate the procedure of taking averages of the soft gluon creation and annihilation operators in 
the valence projectile state. The operator Ui[S,^] corresponds to our operator S^y^ where all the valence charge 
densities and the soft gluon operators are rotated by the matrix S, while the operator Uf[2,^] corresponds to Sly^ As 
discussed in the previous section the operator is responsible for the emission of soft gluons in the incoming wave 
function. These gluons scatter eikonally while propagating through the target. This is precisely the role of Ui[S,^] 
in the formalism of 1431 . On the other hand the rightmost factor fly^ in eq. pil.3p is responsible for the final state 
emissions, which in [43] is achieved through the introduction of Uf[2,^]. The random noise variables ^ reproduce the 
contributions of the averages of the type aa^ , when both the creation and the annihilation operators arise from the 
expansion of fip or fJp in the same amplitude (to the right of the operator O in eq. pil.3p ) or complex conjugate 
amplitude (to the left of the operator O in eq. pil.3|) ) . The variable S reprodices the contraction whereby a from the 
conjugate amplitude is contracted with from the amplitude. These contractions are affected by the presence of the 
operator O, and for the diffractive cross section exclude the contributions of the gluon operators in the gap. 



D. Diffractive observables in the dipole model limit 



In the previous section we have discussed variety of diffractive observables in the JIMWLK/KLWMIJ approach. 
The evolution of these observables with respect to different rapidities in the process can be expressed in terms of a 
functional of two variables Z[S,S]. The evolution of this functional across rapidity gap is given by the sum of two 
independent JIMWLK Hamiltonians, while the evolution across intervals where no restriction on final states is imposed 
is with the mixed JIMWLK Hamiltonian eq. (jIII.17[) . The difference between the various observables is essentially in 
the initial conditions for the evolution. This is a conceptually simple result. However just like for the case of the total 
cross section, the structure of the functional evolution is complicated, and in fact even more complicated, since the 
number of degrees of freedom is doubled. 

It is thus useful to understand these observables in the dipole model limit which puts the evolution in terms of 
simple rather than functional differential equations. Some of the observables have been discussed before but we present 
the discussion here for completeness. As discussed in Section 2, the evolution of the cross section of a projectile is 



17 



very simple. In terms of the cross section at initial rapidity 1 — S^^[s], the cross section at rapidity Y is given by 
1 — I]^'^[sy] with sy satisfying the Kovchegov equation eq. (|II.4ip . For the elastic cross section ea. (|III.13p the dipole 
limit is clearly given by 

nE ((1 - sr[«r])>T - n (ni.53) 
where again sy is the solution of the Kovchegov equation with the same initial condition as for the total cross section. 



1. Projectile elastic - target inclusive 



Starting with eq. (jIII.14p . in the dipole limit we obtain 

Nf = ((1 - ^r[sY]f)T (III.54) 

which for a single dipole projectile becomes also Ty- 

Thus in the dipole limit N§ = Nf . This is obviously the consequence of the target factorization approximation. 
Eq. piI.14p involves target average of the product of two identical dipoles sitting on top of each other in the transverse 
plane. As discussed in Section 2 this is the situation in which we expect the assumption of the independent scattering 
of the two dipoles to be maximally violated. For possible phenomenological application it is therefore worthwhile to 
relax the target mean field approximation by introducing as an independent degree of freedom the correlator of two 
dipoles y)s{u, v))t)- Note that when the two dipoles are identical this is the same as the scattering cross section 

of a gluonic dipole. This correlator probes target field fluctuations as discussed in Refs. [25l. Bsl |49|. fsBl . Following 
proposal of Ref. one can construct a Gaussian distribution for target weight functional. The mean value of the 
Gaussian would determine (s)t while the variance would be given by (s^)t. Within this setup we would have 

A^f =((1 - srM)')T (111.55) 

where the averaging is performed over the initial Gaussian distribution which characterizes the target [5^. The variance 
can be fixed by the data on the elastic scattering at lower energy. 



2. Elastic target scattering inclusive over the projectile 



We refer to ea. pil.isp . Since the averaging over S and S factorizes, it is easy to see that in the large Nc limit the 
"composite" dipole made of 5^5 factorizes into the product of two "elementary" dipoles 

-l-(tr[SHx)Six)SHy)Siy)])T = {^tT[Six)S\y)])T {^tT[Siy)SHx)])T (III.56) 

Since each target weight function evolves according to JIMWLK equation, this means that the "elementary" dipoles 
evolve according to the Kovchegov equation. Thus in the dipole model limit we have 

A^l = 2Ty - 1 + E^^[4] (III.57) 

with Sy satisfying the Kovchegov equation eq. pi.4ip with the initial condition 

sy^o = s. (III.58) 



3. Projectile diffraction with elastic target scattering 

This observable in the dipole limit has been discussed by Kovchegov and Levin [soj. Consider ea. pil.3ip . In the 
dipole limit the evolution of with respect to the diffractive interval at fixed Ygap is given by the simple dipole 
evolution. Moreover at Yp = this observable reduces to the one discussed in the previous subsection. Thus to obtain 
j\fD,Yp total rapidity Y we can start evolution at Yp = 0, evolve over Ygap to get the observable in the previous 
subsection and subsequently evolve it according to the dipole model over Yp. Thus we have 

Tvi?'^- = 2ry - 1 + j:r[st,Y.J (in-59) 
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where Sy^ y obtained by solving the Kovchegov equation with respect to Yp with the initial condition 

4.=o,y,., = 4,., ■ (111.60) 

The derivation in ^3^] is given for the projectile which is a single dipole. Eq. (jIII.59p is the generalization for 
an arbitrary initial projectile wave function E^^[s]. Diffraction via the Kovchegov-Levin equation was extensively 
investigated in Ref. [57| . 



4- Projectile diffraction inclusive over the target 



This observable is similar to the one discussed in the previous subsection. The relation between the two is the 
same as between the totally elastic scattering and projectile elastic, target inclusive cross section. Indeed, examining 
eas. piI.24IIII.25P we see that 



^2Ty-l + (I]^^[.^^_y_])T (III.61) 



Here 



(III.62) 



and the subsequent evolution of Syp y '^^^^ according to the Kovchegov equation. Thus clearly if we assume 

target mean field approximation for the averaging over the target in eq. pil.61|) we return to the Kovchegov-Levin 
observable of the previous subsection. However the mean field approximation for eq. pil.6ip is maximally violated, 
as we have to average products of at least two dipole operators at the same point. The proper way of calculating 
this observable therefore is again to have an ensemble of configurations for the target averaging. The evolution then 
has to be performed for each element of the ensemble. For each element of the ensemble we will obtain the analog of 
Kovchegov-Levin 7V^'^^[s] which then has to be averaged over the target ensemble. In this sense we have 

TVf '^^ = {N^^^^)t (in.63) 
and we expect the target ensemble averaging to give important corrections. 



5. Double diffraction. 



In this section we discuss the fate of the double diffractive cross section in the dipole model limit. We start by 
considering eq. pil.4ip . Our first observation is that for any finite gap only color singlet intermediate states contribute 
in the sum over P' in cq. pil.4ip if the initial state is a color singlet. The physical reason for this is very simple. 
Recall the definition 

E^^^; ^ {P.\n^l^ s n^^JP^) . (in.64) 

Let us suppose that \Py) is a physical color singlet state localized in the impact parameter plane. As discussed in detail 
in [111 the JIMWLK evolution has the property that the wave function of the state spreads in the impact parameter 
plane. However if the state is color singlet this spread is rather mild - the long distance tails that are generated by 
the evolution decrease as ^ . Thus the probability density to find partons in such an evolved state decreases towards 
the periphery as Thus after the evolution the state \Py) is still localized with all the probability concentrated at 
central impact parameters. On the other hand for a color nonsinglet state |P^) the situation is radically different. 
The Coulomb tail generated by the evolution decreases only as jij and the probability density decreases only as 
Thus after any finite evolution interval all the probability for such a state is concentrated at spatial infinity. It thus 
follows immediately that an overlap of an evolved singlet state and an evolved colored state vanishes no matter how 
small the evolution interval is. The presence of the S- matrix operator in eq. pil.64p does not affect this conclusion, 
since the action of S is completely local in the transverse plane. 

It is easy to put this argument into more technical terms. Acting by the JIMWLK Hamiltonian on the S'-matrix 
element of eq. pil.64p we find that the infrared divergences do not cancel in the virtual part. This (negative) divergence 
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exponentiates for a finite evolution interval and the matrix element vanishes. In Appendix C we present this calculation 
for a matrix element between a singlet and an octet dipoles in the dipole model approximation. 

We conclude that as long as Ygap ^ 0, the colored intermediate states have to be omitted from ea. pil.4ip . 

\ singlet P' 

+ J2 J DSWl[S] St ^^'[5] ^P'P[S] 5{Y,,,) . (III.65) 

colored P' 



The form ea. pil.65p is suitable for taking the dipole limit since all the elements in it are singlets, and therefore can 
be taken to depend on the dipole degree of freedom only, E^-^ [S] = [s] . 

Let us consider the evolution of the cross section with respect to the rapidity gap. Since each factor E evolves with 
the dipole Hamiltonian ea. pi.38[) . and since the dipole Hamiltonian is first order in the derivative with respect to s, 
we have 



gap 



f DS wl [s] '^-^ [.] ( 1 - [.] - M + E ^yZ' t-^] ^yZ W ] (ni.66) 

y singlet P' J 



Note that the dipole Hamiltonian is Hermitian with respect to the proper integration measure DS, and therefore we 
can integrate the Hamiltonian by parts and put it on . We thus have 



\ singlet P' 

<1 (111.67) 



We see therefore that evolving the cross section with respect to the rapidity gap is the same as evolving it with respect 
to the target diffractive interval. This establishes that the cross section does not depend separately on Ygap and Y^, 
but only on the sum Ygap + Yt ■ 

The width of the projectile diffractive interval is not essential for the argument. Clearly, we can equally well allow 
the projectile to diffract in any finite rapidity interval Yp. The cross section then depends separately on Yp and 
^gap + Yt- Thus the double diffraction in the dipole model is equal to the single diffraction. It is also clear that 
imposing further gaps and/or diffractive intervals at intermediate rapidity does not change the result. The diffractive 
cross section depends only on two rapidity variables: the diffractive interval of the projectile and the total distance 
in rapidity between the diffractive remnants of the projectile and the valence rapidity of the target. 

To restate our conclusion, we find that to define diffractive scattering within the dipole approximation it is sufficient 
to sum over color singlet states in some rapidity interval Yp on the projectile side and require an arbitrarily small 
gap below this interval. This automatically ensures that there are no gluon emissions on the target side of the gap. 
We note that in [i^ the diffractive scattering was defined indeed simply by summing over color singlet intermediate 
states on the projectile side. 

Note however that the argument does not extend beyond the dipole model, or rather beyond the leading Nc 
approximation. It was crucial for our proof in eq. pil.67|) that we could represent the evolution with respect to the 
width of the gap as the single Hamiltonian acting on the sum of the projectile averaged S - matrices. This holds in 
the large Nc limit, as the Hamiltonian is linear in the functional derivative. However the full JIMWLK Hamiltonian 
is a quadratic functional of the functional derivative with respect to S. Thus the evolution of the last term in the first 
line in eq.(jTTL65| can not be represented as the JIMWLK Hamiltonian acting on the product E^^^ [^JE^ ^[S] even 
if the states P' are color singlets. Consequently the evolution with respect to the gap cannot be traded for evolution 
with respect to Yt- We conclude that the subleading in 1/Nc terms in the JIMWLK Hamiltonian are responsible for 
the double and multiple diffraction processes discussed in the previous section. 

This concludes our discussion of diffractive processes. We now turn to examples of other types of observables. 
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IV. SCATTERING WITH MOMENTUM TRANSFER 



In this section we consider scattering processes with fixed transverse momentum transfer. All of the observables 
considered here are inclusive with respect to the final states of the target. 

A wave function of a probe with definite transverse momentum p can be written as: 

|Pp) = ^rp(r) = J d^e'P^'i'ix) (IV.l) 

where b = ^i/n is the impact parameter of the configuration of n gluons at transverse positions x^, and denote 
the relative distances between the gluons. 



A. Elastic scattering with momentum transfer 



We take the initial projectile state to have transverse momentum zero and the out state to be the same state but 
with transverse momentum q. The operator observable that is being measured has the following form: 

= \P,) (P,\ = [ <fh<fh 6^(^)0-^) in) . (IV.2) 



(2 7r)2 ' ^ «' (2 7r)2 

This corresponds to putting the separation scale Yq = Y , that is the momentum transfer is fixed for the valence part 
of the projectile wave function. For the expectation value of the observable we obtain 



(IV.3) 



N""-' = ^ j DSW^[S] J ShSve^" (l - ^l''''[S{x)]) (1 - + v)]) . 

The Y evolution of N^-'^ is given through the evolution of with the Hamiltonian jjJIMWLK _ 
In the dipole limit, assuming the target mean field approximation we have 

with 

N{q,Y) - / d% (1 - ^^''[syM) e^'"- (IV.5) 



Since the dipoles in the final state are displaced with respect to the dipoles in the initial state, the mean field 
approximation is not suspect in this case. 



B. Total cross section with momentum transfer 



We now sum inclusively over the final states of the projectile with momentum transfer q. 

N''" = / DSW^S] I dHd^e^^" Y: ('5^^' - ^l^'^'iSix)]) {s^^' - J:^'^[S{x + v)]) (IV.6) 

where the sum over P' runs over all final states of the projectile. The Y evolution of A'^^^'^ is again given by the 
evolution of W"^ with H^^^^^^^ . 

Let us now consider this variable in the large Nc limit. It is easy to see that even in the large Nc limit and 
even assuming that the incoming state P contains only singlet dipoles, the observable ea. pV.6[) can not be calculated 
without introducing quadrupoles. Consider for example a non-forward scattering of a single quark dipole. There are 
two intermediate states that contribute to the scattering, the color singlet and the color octet. For the color singlet 
state as usual we have 

^^""{x^y) ^^tr[Sp{x)Sl{y)] (IV.7) 
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For the (normalized) color octet intermediate state P' we have (see Appendix C) 



^0 



■tr[SF{x) r"4(y)] 



Summing over the singlet and octet states with equal weights we find 

(fhcPv 



a{r,q-Y) = / DSW^iS] 



'tr 



a{r,q-Y) = 2TY5\q) - 



{2n] 



(l - SF{x)Sl{y)'^ (l - SF{x + v)SUy + v) 
2 , (fbd^ve'-'"'' (1 - q{x,y,x + v,y + v)) 



with 



q{x,y,x + v,y + v) = Sf{x) Sy{y) Spiy + v) Sf{x + vf 



(IV.. 



(IV.9) 



(IV. 10) 



(iv.ii) 



The last term involves the quadrupole scattering probability and is not suppressed by powers of relative to the 
first term. In fact it is easy to see that for an arbitrary projectile made of dipoles only we have 



N'^^ = 2TY5\q) - J DSW^iS] J d^d^ve'-^l - J:r[q{x,y,x + v,y + v)]) 



(IV. 12) 



where T is the total cross section. Thus we see that in the large Nq limit it is not sufficient to specify the average 
of the dipole amplitude in the target wave function, but one also needs to specify the quadrupole amplitude. The 
JIMWLK evolution of in pV.12|) can be integrated onto As discussed in Section 2 the evolution of T,^^ is 

then 



(IV. 13) 



where qy is the solution of the quadrupole evolution equation pi.52[) derived in Appendix A with the initial condition 

qY=o = q- 



C. Total momentum transfer within a fixed rapidity interval. 

In defining this observable we require that the total momentum q is transferred inside the rapidity interval Yp on 
the projectile side. The generalization of eq. (|IV.6p is straightforward 

iV^^^-^ - / DSW^_y^[S] J ^^e-^ ^ (S^^' J:1^/[S{x)]) {S^r>' _ ^p;p^s{x + V)]) 



2TYS^{q) 



DSW^^yA^] / d%d^ve''"i (1 - ^^^[S\x)S{x + v)]) 



(IV. 14) 



This variable can be evolved both with respect to the total rapidity Y and the rapidity interval Yp. A convenient 
representation for this purpose is 



j^IXp,q ^ 2TY5^{q) 



(2 



^ j DSDS j dHd^ve'"^ Zy.Yp[S,S] (1 - i:^P[S\x)S{x)]) . (IV.15) 



The weight functional Zy,Yp is found by solving 

d 



dYp 



Zy.Yp [S., S] — — H^IS] S] Zy,Yp[S, S] 



with the partial derivative taken at fixed Y — Yp and the initial condition for the evolution 

Zy-Yp.Yp=o[S,S\ = W^_Yp[S] 6[S{x)-S[x-v)] 
with Wy^Y ["^l evolving according to the JIMWLK equation. 



(IV. 16) 



(IV. 17) 
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D. Diffraction with momentum transfer 



The simplest observable of this type is elastic projectile scattering with total momentum transfer q and rapidity 
gap Yp. For this observable we have 

N"^^'^"^- - J DSW?_y^JS] J ^e'- (l - 4f„f[^(x)]) (l - + -)]) (IV.18) 

The evolution of this observable with respect to total rapidity Y at fixed Ygap is still given by the JIMWLK evolution 
of W^^. 

Finally we can ask for total momentum transfer in a diffractive process where the projectile diffracts into the 
rapidity interval Yp. Combining the results of Section 3 with the earlier discussion in this section we can write 

= 2Ty5\q)- j DSW^[S] J ^^e--^ (^1 J] E^f/^ ^^^.^^iYp; S{x + v)]^ (IV.19) 

and the summation has exactly the same meaning as for the diffractive process discussed in Section 3.3.2 and 3.3.3. 
Following the same logic as in Section 3.3.3 we can rewrite it as 

= 2TY6\q) - ^ I DSDSZy,^y,JS,S] J <fhd\e-"'^ (l - ^pP[S\x)S{x)]) . (IV.20) 

s before is 

Zy,y,^, [5, 5] - - H^[S: S] Zy.y,^, [S, S] (IV.21) 

yp=o,y,„, [S, S]^ -H2 [S; S] Zy,=o,y,„, [S, S] . (IV.22) 



(2^) 

And the evolution of the functional Z as before is 

d 
dY^ 

and with respect to the rapidity gap as 

d ^ 
BY ' 

^ ^ gap 

Like in Section 3.3.3 the initial condition for the evolution with respect to Ygap is 

Zy,=o;Y,^,=o[S, S] = W^[S] 6[S{x) - S{x - v)] (IV.23) 

This have to be evolved first with respect to the gap to Ygap and subsequently with respect to the diffractive interval 
to Yp. 

To get the large Nq limit for this observable we have to understand the evolution from the point of view of 
rather than Z. At zero gap we simply have the observable of the previous subsection and we need to know the form 
of Sy^. This is clearly obtained by 

Ei^f = (IV.24) 
with s'^' solving the Kovchegov equation with the initial condition 

^Yp^oi^' y) = y,x + v,y + v). (IV.25) 

Note the difference between Eq. (jIV.13|) and Eq. pV.24p . In Eq. (jIV.13|) the evolution is that of the quadrupole while 
in Eq. pV.24[) the dipole evolution with the quadrupole entering as initial condition only. The subsequent evolution 
across the gap is with two independent JIMWLK Hamiltonians with respect to S and S. We also know that the octet 
states do not make it across the gap in the large limit. This means that for the purpose of the Ygap evolution we 
can write 

q{x,y,x + v,y + v) = s{x,y) s{x + v,y + v) (IV. 26) 

and evolve each s according to the Kovchegov equation. All said and done, the dipole model observable is obtained 
as the target average of 

TVW- = 2TY5'{q) - J d'bd\{l-J:r[4'p,Y,J) (IV.27) 

with Y {x, y) evolved with Kovchegov equation with respect to Yp from the initial condition Sy^^g y (x, y) = 
SYgap{xTy)sYg^p{x + v^y + v) with SYgapi^iV) evolved by the Kovchegov equation from the initial condition 
SYg^p=o{x,y) = s{x,y). 

This concludes our discussion of transverse momentum transfer. 
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V. INCLUSIVE GLUON PRODUCTION 



The last observable we consider is the inclusive gluon production. Within the dipole model this has been discussed in 
[40| . In [i^ the inclusive gluon production was calculated without the dipole approximation, but the rap idity evolution 
although implied was not exphcit. Single gluon production was also discussed in Refs. [ssl [sol [60l [6ll |6^. [63j 

We are interested in a differential cross section du/dydk'^ for production of gluon at rapidity Yq and transverse 
momentum k. At this rapidity the observable is given by the expression [42]: 

Og = A(ojf]^(l - S)nl n{k,y) ny{l - S^)^l\Qa)\y=Yo (v.i) 
d y 



^e" A 

n{k,y) = / dk~ a\''{k,k-)a1{k,k-) (V.2) 

The target W'^ is evolved with the Hamiltonian ff-'iMWLK ^^le rapidity Yq, while the projectile is evolved 
with H^LWMij rapidity Y - Yq 



da 



J DS Wl [S] J Dpp W^_y^ [pp] dg{Yo) (V.3) 



dYo dk^ 

For KLWMIJ/JIMWLK evolution the operator fly is known explicitly ea. pi.3[) . and Og can be computed explicitly 

0,^f fllle^H^-^) Q1{z)Q'i{-z) (V.4) 

J Z TT Z TT 

For a complete evolution operator which includes Pomeron loops the observable will depend also on the multigluon 
production amplitudes Qn and it is not known at present. 

We now present a short derivation of Eq. (|V.4p . From Eq. pi.Sp we have 

at{k,k-)Cy = Cy [at{k,k-) - -^K[PP]) (V.5) 

and 

[a^{k,k-) - -^K[PP]) (1 - S^)Cl\Qa) = {a^{k,k-) - -^K[pp]) Ct|0,) 

-S^ (^S^'atik,k-) ~ -±=bnSpp]^ Cl\Oa) = -^S^C'liS'^Htipp] - bnSpp])\Oa) (V.6) 

To leading order in the coupling constant the operator C commutes with b. Thus we obtain 

Oy = [S'^'blipp^k] - bnSpp;k]) {S^H'^[pp-k] - bnSpp;k]) (V.7) 
This coincides with Eq. (|V.4p if equivalently written in coordinate representation 

Oy^ f e^'^^^-'US'"'bt[pp;z]^bnSpp:z]){S^^bnpp;z]~bt[Spp;z]) (V.8) 

If we keep the nonlinear terms in the expression for the "classical field" b{p), Eq. (jV.7[) and Eq. (jV.Sp provide a 
generalization of the results of [l^] to include some non- linear effects in the projectile wavefunction (Pomeron loops). 
In that case for consistency we have to use the JIMWLK+/KLWMIJ+ evolution rather than JIMWLK/KLWMI J 
discussed in the bulk of this paper. 

Back in the JIMWLK/KLWMIJ limit the operator Og is: 



O - ^ 



_ j pa p6^(y) r^^ _ s.sl- si + sl\ (V.9) 



Jz,z Jx.y {Z - X)^ {z - yf 
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As has been noted before [l^], [1^ the operator Og Eq. (jV.9p respects the projectile-target factorization. For 
evolved with the KLWMIJ Hamiltonian the correlator {pp{x) p''p{y))p satisfies the BFKL equation. Thus for color 
singlet projectile and target 



da 



dYn dk^ 



{z - x)^ {z ~ y), 
[z — xY [z — 



with {Tx^y)T = N'2_i (tr[»S'(x)5''l^(^/)])T standing for a gluonic dipole scattering amplitude and n defined as 



(V.IO) 



{x,y;Y-Yo) = / DppW^_y^^[pp] pU^) p'^piy) . 



(V.ll) 



This can be written as 



da 



dYo dk^ 



P-.q 1^ 



{k+p)z + {q—k)z —px — qy 



X nP{x,y;Y-Yo) m,y)Y„ + {T,,z)y„ - {T,,z)y„ - {T^.,yh 



(V.12) 



This expression does not assume the dipole model limit nor is it restricted to an initial state consisting of a single 
dipole. The wave function of the projectile enters only through the initial conditions on the evolution of . 

In terms of the Fourier transforms 



- iki X - ik2y P 



n {x,y;Y - Yq) ; 7^,(^1, fc2) = 



- i ki X — i k2 y 



{Tx. 



y)Yo 



Eq. (rvTT2)) takes the form 



da 



dYo dk^ 



TT 



K2^2{k; q, p) nY_Yo{Pi<l)TYo{-k - p, k - q) 



(V.13) 



(V.14) 



with the vertex 



K2^2{k; q, p) 



qi ^ ki 




ki Pi 


q2 ^2 




jj2 p2 



(V.15) 



If either projectile or target is assumed to have translational invariance in the transverse plane, the constraint p = —q 
is authomatically imposed. In this approximation the result reduces to the standard kt factorizcd form. 

For future applications we note that it is sometimes useful to perform the integral over p^ in Eq. (jV.Sp . This 
averaging procedure as always turns into T,[S] and any additional factor of p^ present in Og into right or left 

SU(N) rotation generators Eq. (jll.32|l . To this end it is useful to temporarily set the second factor of S in eq. (jV.3|l to 
S. Some algebra then gives 



da 



dYn dk^ 



Jk(z - z) 



J DSDSS[S~S]W^^[S] Q^{S,z)QfiS,z) I^y-YoIS^S] 



(V.16) 



This representation is similar to the one used in Sections 3, 4 and 5 for other semi inclusive observables. 
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APPENDIX A: APPENDIX: EVOLUTION OF THE QUADRUPOLE OPERATOR 



In this Appendix we consider an evolution of a quadrupole operator [4ll . |47 

Qx,y,u,v — tl"{Sx Sy Su Sy] 



(A.l) 



The evolution of q follows from the action of the JIMWLK Hamiltonian on q (second equation of the Balitsky's 
hierarchy [12]). Let us introduce three kernels. The Weiszacker- Williams kernel 



K. 



x,y,z 



{x - z), {y - z), 
{x - z)2 (y - zf 



The dipole kernel 



The "4 ^ 4" kernel \2A U 



L 



x,y,u,v;z 



{x - z)i 



iy 



K. 



{x - yf 



[x — z)2 {y — zY 



^'"'^ {x - zY {y - zf 

ju - z)i ^ {v - Z)i 
(u — z)^ (v — z)^ 



(A.2) 



(A.3) 



(A.4) 



To derive the evolution of q we follow the same procedure as for the derivation of Kovchegov equation. Namely we 
take the evolution equation for four Wilson lines and factor its right hand side using the large Nc factorization. The 
result is [13 



dq: 



x,y,u,v 

dY 



— 71 / \-^'^x,y\z -^^u.v.z ~^ Lx^u,v,y\z\ Qx,y.u,v -^x ^y.u^v.z ^x,v '^y,u ^x.v .u^y.z '^x.y '^u,v 

~^ ^x.v,u,v\z Qx,y,u,z ^z^v ~^ ^x .y ,x ,v\z Qz,y,u,v ^x,z ~^ ^x .y ,u,y;z Qx,z,u,v ^z.y H~ ^u,y ,u,v;z Q_x.y,z,v ^u,z} 



(A.5) 



Note that as opposed to the Kovchegov equation, Eq. (|A.5[) is linear in q. It is however coupled to s whose evolution 
is nonlinear [41j. In order to compute the energy behavior of the quadrupole, one first needs to solve the Kovchegov 
equation for dipoles, and then solve the linear equation for the quadrupole coupled to dipoles. Note that due to the 
inhomogeneous term in Eq. (|A.5|) even for the initial condition g = a nonvanishing quadrupole is generated by the 
evolution. 

We now comment on the physical meaning of the quadrupole operator q. Consider a projectile which consists of 
two dipoles at points x, y and u, v. The normalized wave function of such a projectile (in terms of the creation and 
annihilation operators of the " quarks" and " antiquarks" ) is 



\{xy),{uv)) = ■^^qt{x)q,{y)^qj{u)qj{v)\Q); 



(A.6) 



where i, j are color indices in the fundamental representation. After propagating through the target fields the state 
emerges with the wave function 



\ixy),{uv))out = ^ J2^S{x)qix)],[SHy)q{y)]^ Y.lS{u)qiu)], [SHv)qiv)], \0) 

i 3 

It is easy to calculate the overlap of this outgoing wave function with the incoming one. 

{{xy),{uv)\{xy),{uv))out = s{xy) s{uv) . 



(A.7) 



(A., 



One can also calculate the overlap of the outgoing state with the two dipole state where the quarks exchanged their 
antiquark partners 



{{xv),{uy)\{xy),{uv))out = q{xyuv) 



(A.9) 
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Thus the quadrupole operator is the antiquark exchange amphtude. For the single dipole pair of ea. (|A.6|l the 
quadrupole does not contribute to the total cross section, or forward scattering amplitude. However a general color 
singlet state with two quarks and two antiquarks is a superposition of two possible dipole pairs 

\singlet, xyuv) = a\{xy), (uv)) + f3\(xv), (uy)) (A. 10) 

The forward scattering amplitude for such state is 

S = aa* s{xy)s{uv) + /3(3* s{xv)s{uy) + — — [a* (3q{xyuv) + (3*aq{xvuy)] (A. 11) 

Nc 

Thus the quadrupole contributes 1/Nc correction to the total cross section. In the leading order in 1/Nc therefore the 
quadrupole contribution is absent. However for a state with n dipoles there are dipole pairs which can exchange 
antiquarks. Thus the quadrupole contribution can become important (depending on the exact wave function) already 

1/2 

when the number of dipoles n ~ NJ . On this note we also mention that the equation for q contains an inhomogcneous 
term. Thus even if q vanishes at initial rapidity it is generated through the evolution. This in particular means that 
the dipole model limit is not recovered from eq. (jH.47p by setting g = but rather by dropping the q - dependence in 
the function S[s,(7]. 



APPENDIX B: APPENDIX: ATTRIBUTING FINAL STATES TO THE TARGET HILBERT SPACE 



In the body of this paper we have avoided introducing intermediate states in the target Hilbert space. Nevertheless 
this can be done. In this way we can relate the functional Z[S, S] introduced in Section 3 to the density matrix of 
the target. 

We take the target wave function to be boosted to rapidity Y. To discuss the resolution of identity on the target 
Hilbert space it is convenient to introduce the basis of eigenstates of the operators Sy = Pexp{i j_ dx~T°-a°'{x, x")} 
for allO<y<Y: °° 

Sy \{Sy}) = Sy \{Sy}) t^^'lS] \{Sy}) = \{Sy}) (B.l) 

Any target state (Ty| which depends on all intermediate rapidities < y < Y can be expended in this basis: 

\Ty) = J DSy Q^[Sy] \{Sy}) J [S y] ^ [S y] 0*^" [S y] J {Ty \ Ty) = j DSy [S y] = 1 . (B.2) 

In this equation the functional integral is over Sy for all < y < 1". This is defined simply as the integral for all x~ 
dependent a{x~): DSy = Tix-Da'^{x~). 

Just like for the projectile off diagonal matrix elements it is easy to see that both Wy as well as the target off- 
diagonal matrix element Oy evolve with the Hamiltonian jjJIMWLK ^ This is a consequence of Lorentz invariance 
for any element of the ^-matrix. Resolving identity on the target side in this basis for an arbitrary observable O 
Eq. plL4)) we get: 

{0)y = J DSyDSyY,[6^>^' - J:l%Py]) (<5^^" - ^^J^,]) 

p'p" 

X el[Sy]e*y^[Sy]OP'P"[Sy,Sy]. (B.3) 

In particular consider an observables which is completely inclusive over the projectile degrees of freedom. This 
includes projectile diffractive observables when we choose the separation scale Yq such that it includes the rapidity 
gap on the target side. 

OP'p"[S,S] ^ 0[S,S]6p,p. (B.4) 
Eq. (jB.Sp can be rewritten in the following form 

{0)y = J DSyDSyJ2{S^'' ^VXo^Sy]) ^S^' ^ - ^^^^y^lSy]) y] 0^.^ [S y] 0[S y , S y] . (B.5) 

P' 
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Taking into account that 



we can write 



Y,^'^y%^[Sy]^i^_y^[Sy\ = S^^^J^t^,] (B.6) 



{0)y = j DSyDSy (l - - s;^^J5t] + S^^^J5t5,]) Z^J^,,^,] (B.7) 

where we have introduced 

zl [5, s] ^ el [s] e*yj[s] o[s, s] (b.s) 

This is in the form Eq. (|III.32|1 . Note that the functional Z includes information both on the target and observable. 
For example for elastic target scattering 

o[s,s] = e*yl[s]el[s] (b.9) 

and we recover Eq. (|III.3ip . With respect to Yq this observable obviously evolves with H'"^-^^^^[S]+ H''^^^^^^[S]. 
On the other hand the evolution of any Yq independent observable follows from the known evolution of the projectile: 

^z^[s,s] = - ^ {Qns] + omY z^[s,s]. (B.io) 

This is precisely the evolution we have found in Section 3. This illustrates that we can put the separation scale Yq 
either above or below the rapidity gap with the same results for the evolution, as expected. 



APPENDIX C: APPENDIX: INFRARED DIVERGENCE OF THE COLOR OCTET EVOLUTION 

In this appendix we derive the evolution equation for the scattering amplitude of a dipole into a color octet final 
state. 

The normalized singlet and octet states defined via "quark" and "antiquark" creation operators are 

\P) = \singlet) = ^ g-.(y) |0) ; \P') = \octet) = ^ q{y) t'' q{x) \G) (C.l) 

The ^'-matrix element reads 

^^-y' ^ \[^tr[Sl{y)T^^SF{x)] (C.2) 

We act with the JIMWLK Hamiltonian on S^J^ . Evaluating separately each of three terms in the Hamiltonian we 
obtain 

/ K^,v,.J'k{u)Ji{v) Sf^^' ^Cf f M,,y,,E^f (C.3) 

'J u,v,z J z 

I K^,v,. Jl{u) Jl{v) ^Cf I [K,^,^. + K,,,,,] + ^ / [K,,v,^ + Ky^.A Ky' (C.4) 

J u,v,z J z ^ -^^ J z 

— N I ^Kx^x.z ~ Kx,y,z\^zy ^x,z + [Ky_y_z ~ -^y,x,z]^xx ^z,y^ (C-^) 
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Cf = [N'^ — 1)/2N. As opposed to the equation for the singlet transition ampUtude, there is no cancellation of the 
infrared divergencies as z —^ oo between the virtual terms (J^ Jn and J^) and real term: 

+ N ^ {[X.,.,. - K,,y,^] Sf/' S,^, + [Ky^y^, - Ky^,^,] S , ^y^ (C.6) 

Consequently the total action of the JIMWLK Hamiltonian on is divergent and negative as the virtual terms 

enter the evolution equation with the negative sign. We thus conclude that '^^y' vanishes after evolution over 
arbitrarily small rapidity interval. 
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